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Methods for calculating the power in the quantizing noise on digital 
transmission facilities have been knoiim for some time. A more difficult 
but unavoidable problem is the effect that this noise has on data signals 
intended for analog transmission. This paper demonstrates that to assume 
that the noise will behave as a white Gaussian noise process will always 
(except for a simple factor) yield an upper bound on the probability of 
error when no companding is present. We assume that linear detection will 
be used, as for a PAM system, and the result is true whether or not filtering 
or demodulation is involved. Results are illustrated by applying them 
to a model of an existing VSB modem luhereby the additional degradations 
resulting from data set imperfections are included as added baseband noise. 

A modem operating perfectly would make no errors at all at the higher 
transmission levels. For example, with no companding, a set with an eight- 
level eye closed by even 80 percent would not yield errors for input powers 
down to — 15 dBm. Thus quantizing noise is not a basic limiting factor in the 
error rate for all input levels. A similar rigorous theory is not available for 
compandored systems, but for special situations reasonable estimates can 
be made. For logarithmic companding and eight-level VSB transmission, 
worst case estimates indicate error rates about 10~ 6 for one link of Tl 
carrier. 

I. INTRODUCTION AND SUMMARY 

The Tl carrier system is a digital transmission scheme for analog 
signals. 1 Even though the digits in the coded bit stream might be 
transmitted without error, when the analog signal (which may in 
fact be a data signal designed for analog facilities) is reconstructed 
at the receiving terminal, quantizing noise is inevitably added and 
can be large enough to cause errors in the customer's data. 
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We show that, under some simple constraints between sampling 
rates and band widths which are satisfied in practice, and independent 
of the particular data signal used, an upper bound on error rates is 
obtained if the quantizing noise is assumed to be a white Gaussian 
process of power A 2 /12 and bandwidth 1/27 1 , .* We assume that linear 
detection will be used, and the result is true even if additional filter- 
ing is done (as one might do with a receiving filter). And it is true 
whether or not a demodulation process takes place. Using the model 
in Fig. 1 for the digital transmission scheme, results are specialized 
to obtain error rates for eight-level VSB transmission (Fig. 2). 
Imperfections of the data set are included as added baseband noise. 
If it were not for these imperfections, error free transmission would 
result over an appreciable range of power levels (see Table I). For 
a logarithmic compandor and VSB data, even using worst case esti- 
mates, the error rate for one link is quite low, about 10~ fl . 

II. QUANTIZED TRANSMISSION SCHEME 

Let us consider a transmission scheme for a single channel that, 
for our purposes, typifies the Tl carrier system. As suggested in Fig. 
1, the signal to be transmitted is assumed not to have any power 
beyond B Hz. The signal is sampled at the Nyquist rate T^ = 1/(2J5) 
and these samples are passed through an instantaneous nonlinear 
device with characteristic v ont = F(v in ). The compressed samples are 
then quantized by a uniform quantizer of step size A, and coded into 
binary sequences. The binary sequences are assumed to be trans- 
mitted without error and the process is reversed: sequences are de- 
coded into pulses, expanded according to the inverse function F^ix) 
and the impulses are used to excite an ideal filter of bandwidth B 
and amplitude Ti .t A receiving filter G generally follows the ideal 
filter and we include this in our description, although it would not 
be part of a Tl transmission system. If the bandwidth of G is en- 
tirely contained in B then one may consider the impulses to excite 
T X G directly. 

To be more specific, we are concerned with two particular com- 
pandor characteristics F(x). One is F(x) = x, that is, quantizing 



* Here A is the quantizer step size and 1/1\ is the sampling rate. Also this 
statement is true only modulo a simple factor given in the text. 

t The amplitude gain of the ideal output filter for the carrier system is chosen 
to be Ti in order that the signal component will undergo no gain relative to 
its sampled values at the transmitter. 
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Fig. 1 — Essential elements of quantized transmission scheme. 
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Fig. 2 — Model of back-to-back transmission of VSB modem. The noise added 
at baseband represents imperfections of the modem. 
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Table I — quantizing noise* 



Mean 
square 
P(dBm) 


n m al 


- 5 
-10 
-15 


0.221 
0.394 
0.698 


-20 
-25 
-30 


1.245 

2.21 

3.94 


-35 
-40 
-45 


6.98 
12.45 
22.1 



* Peak values of quantizing noise as a function of input power P. The noise scale is 
such that a perfect receiver would make no errors for nm., < 1. The power scale is 
such that the quantizer overloads at an instantaneous power of +6 dBm. 

without compandoring. The other case is (in normalized units) 



F(x) = -F(-x) 
HX) ln(l + M ) 



Og^l 



0) 



= 1 x > 1, 

where p., the degree of compandoring, is large. Typically, /x = 100 
for a good approximation to existing devices. 

Finally, when specific values are required, we assume 7 bit coding 
to be used for the quantized samples and use A = 1/63. 

We hasten to add that quantizing noise is not the only degrading 
factor for the existing T x facilities. Apparently mismatch and mis- 
tracking of compressor and expandor cause nonlinearities which are 
responsible for peculiar behaviors of error rate versus signal power 
curves. 2 



in. general theory 

Let us represent the signal l{t) which is to be sampled and quan- 
tized by 

l(t) = x(t) cos u c t — y(t) sin u e t, (2) 

and the sampling wave as 

2 t(t - kT t - r), (3) 
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where the random timing phase is uniformly distributed over the interval 
^ t g T, . The pulse trains representing (2) immediately after 
sampling, compression, and quantization are given by expressions (4), 
(5), and (6), respectively. 

2Z(ftr, + T )8(t -kT x - r); (4) 

2l con > t ,(kT l + r)8(t - kT x - r); (5) 

2Ln, P (fc7\ + r)8(t - kT x - r), (6) 

where 

WfcT. + t) = F[l(kT x + r)] 

is the compressed sample value and Lm P (0 is the particular one of the 
(2 — 1) levels that the quantizer output gives as the value for l aomv (t). 
If we let the subscript "exp" stand for the result of operation of the 
expandor at the receiving terminal, then the impulse associated with 
time (kT t + t) has area 

[Lm^kTr + t)]«„ = IQeTj + t) + e(kT x + t). (7) 

Because the expandor has as its input an estimate of the compressed 
pulse area, the error term t{kT x + t) is not zero but may take any 
value in an interval, that is, 

The spread A(t) that the quantizing error may take is not neces- 
sarily equal to the quantizer step size A when companding is present, 
but is given by the formula (see Appendix A) 

m = TFMT (9) 

In (9), F'[l(t)] is the derivative of the compressor characteristic 
evaluated at that input amplitude of the signal at the time of the 
sampling. The error signal generated at the output of the receiving 
filter is obtained by convolving the impulse train 

MkT i + r)8(t - kT x - t) (10) 

with the impulse response T x g(t) of the receiving filter.* Denoting 

♦Again, g(t) ia associated with the receiving filter of the data set and the 
constant T, ia the gain of the ideal output filter of the carrier system. 
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this noise by n t (t) we have 

n t (t) = TtMkTi + T)g(t - fciP, - r). (11) 

To proceed further we make the assumption that the quantizing 
error «(0 of the output sample in (10) is uniformly distributed over 
the interval 



r-A(o ( m\ t 



thus having mean zero and variance A 2 (t)/12, and that different sam- 
ple errors are independent. Notice that the latter assumption is not 
the same as assuming that different sample values are independent. 

IV. SUMS OF UNIFORM VARIATES 

As (11) illustrates, a basic problem which must be dealt with is 
the probability distribution of sums of independent and uniformly 
distributed random variables. We will obtain an upper bound on the 
tail probabilities of interest by applying the technique of the Chernoff 
bound. 3 - 4 This bounding technique states that if a probability den- 
sity p{x) has a moment generating function (wig/) M(s), where 

[exp (sx)]p(x) dx, (12) 

•CO 

then 

Q = Prob [x Z a] ^ M(s) exp (-sa) r s ^ 0. (13) 

Thus to obtain an upper bound one simply multiplies the moment 
generating function by an exponential, both evaluated at an arbi- 
trary positive s. Actually it is known that there is a best s to choose, 
and it is that one, if it exists, which satisfies the equation 

4- In M(s) = a. (14) 

ds 

Equation (14) assures a stationary value for the right side of (13) 
and it can be shown that such an s in fact minimizes M(s)e~'°. 

For example, for a Gaussian variate of mean m and variance a 2 , 
the moment generating function is well known to be given by 

M(s) = exp[ms + ^J- (15) 
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Thus the best s to choose is, using (14) , 

• = (JL =^- (16) 

Notice that only if a ^ m is this s ^ 0. Thus, as long as a is greater 
than the mean, we have for the Gaussian case 



Q ^ cxp I - — 



(17) 



2C7- 

where (17) results from using (15) and (16) in (13). For the Gaus- 
sian variate under discussion the exact answer is also well known to 
be given by 

where erfc x is the coerror function. 5 In addition, equation (7.1.13) 
of Ref. 5 states that 

exp (-x 2 ) ^ (t)*[x + (x 2 + 2)*][| erfc x], (19) 

and hence the difference between the Chernoff answer (17) and the 
exact answer (18) for the Gaussian case is no more than the multiplicative 
factor (*)*[(p)* + (p + 2)*] where (p)* = (a - m)/[(2)V|. 

We modify this procedure for our problem with the following obvious 
lemma. 

Lemma 1: If G(s) is an upper bound jor the moment generating junction, 
that is, M(s) ^ G(s) for all s, then 

Q ^ e- a G(s), s ^ 0. (20) 

In particular, a positive s = s which satisfies 



y- In G(s) 
ds 



(21) 



is legitimate. 

Next consider a random variable x which is uniformly distributed 
over [-A/2, A/2]. The variance of this variable is A 2 /12, and it has a 
moment generating function M{s) unlt 

. , sA 
smh — „ / * 2n 

• 1/(s) --s*--£(f) ^Tljl- ^ 

9 
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Now the nth term of the sum in (22) is positive and upper bounded by 



YlTl. 

/ 6j n! 



Hence 



M(s) unif :g 2 



sAY 1 
2/ G 



1 



n\ 



- } = exp 



12 2 



(23) 



Thus we have shown that the moment generating function of a zero- 
mean uniform density is upper bounded by that of a zero-mean Gaussian 
having the same variance.* If the uniform variable has mean m the 
theorem is still true if we use instead the moment generating function 
of a Gaussian with mean m. 

We are now ready to write down a whole class of random variables 
which have moment generating functions upper bounded by those of a 
Gaussian of the same variance. Suppose the result is true for two inde- 
pendent random variables, x and y, of variances &j. and o\ , namely 



M T (s) ^ exp 



M v (s) ^ exp 



m 

[¥]■ 



(24) 



Then using the theorem that the moment generating function of a 
sum of two independent random variables is the product of their 
individual moment generating functions, we have 



M x+il (s) = M x (s)M v (s) ^ exp 



= exp - (a-l + al) 



scr x 



exp 



*[f] 



where 



2 , 2 

n, = o-* + 0"* 



is the variance of (x + y) . Thus the moment generating function of 
a sum of any number of independent uniforms of arbitrary means 
and variances is upper bounded by the appropriate Gaussian one 
(same mean and variance as the sum), and thus use of (17) through 
(21) provides a rigorous upper bound for tail probabilities of the sum. 



* A similar theorem was discussed by, Saltzberg for the case of equally spaced 
delta functions. We have followed his method of proof here. 
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V. FURTHER ANALYSIS WITHOUT COMPANDING 

When no companding is present, the independent random variables 
e(&7\ + t) have variance A 2 / 12 and the variance of the noise (11) is 

<r\ s < n «(/)> = T\-^> jE g\t - M\ - r). (25) 

We evaluate the infinite sum in (25) by using the Poisson sum for- 
mula, namely 

t t t& - *T, - r) - i _£ exp [ **"%" T) ](?,(« I) , (26) 

where G 2 (w) is the Fourier transform of g 2 (t). Now since the band- 
width of the filter G is assumed not to exceed 1/27 1 ! Hz, only the 
m = term of (26) contributes and we obtain 

a] = T.-^-GM. (27) 

Equation (27) implies that the noise power measured before the 
receiving filter is A 2 /12. This result has been obtained by Bennett 7 
who also showed that the spectrum of this noise is flat across the 
band. Further, equation (27) is consistent with filtering white noise 
since 

G 2 (0) = f g\t) dt = ^f\ G(oi) | 2 dec. (28) 

An important fact about (27) is that the received passband noise 
power without companding is independent of many properties of the 
signal. Thus it is independent of signal power and multilevel struc- 
ture. It is not independent of rate, however, since this enters implicitly 
into the factor G 2 (0), and likewise it is not independent of roll-off. 
By halving the speed and doubling the number of levels, one decreases 
the noise by 3 dB, but loses 6 dB in noise margin, thus leaving one 
with a net loss of 3 dB in noise margin. Thus it is best to use as few 
levels as possible consistent with given speed objectives, at least if 
the quantizing noise behaves anything like Gaussian noise. 

Let us discuss further some statistical aspects of the quantizing 
noise at baseband. The "line" signal must be demodulated as in VSB 
transmission by multipling the (filtered) received signal by cos u c t 
and eliminating double frequency components. We represent the 
impulse response g(t) of the passband receiving filter G by 

9(f) - 0,(0 cos w c t - g„(t) sin w c t. (29) 
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Further we specialize to the practical constraints* u e > B, and 2tt/T\ 
> 2(w c + Bi). The demodulated noise is 

n,,(t) = I 1 MkT x + r){g x (t - kT t - r) cos fofeT, + co c r) 

+ g v (t - kT x - t) sin [uJcT, + ta e r] } . (30) 

The general expression (29) can be simplified for a VSB receiving 
filter which is symmetric about midband frequency on, and linear 
phase characteristic, by writing 

0(0 = 0i (0 cos co, ^ , 31 v 

oi c — coi = JT/2T, 

where l/T is the symbol repetition frequency. Of course a #* and a f7„ 
may be immediately written down from (31). 

From (30) we derive in Appendix B, equation (32) for the baseband 
variance <rl(t): 

al(0 = \ (| ^)T X [G,M + G,„(0)], (32) 

where G 2l {u) is the Fourier transform of g]{t). We now will show that 
this result is identical to the baseband noise power that would appear 
if flat Gaussian noise of power A 2 /12 were on the line. We do not regard 
this as obvious; in fact it is not true that the signal power at baseband 
is the same as if one had Gaussian noise of the same power and spectrum 
on the line that the signal has. The proof depends on a few simple 
observations. If passband Gaussian noise is represented by 

n(t) = n;(t) cos u c t — n v (J) sin <a e t, (33) 

then 

al = (n 2 ) = (nl) = (nl), (34) 

and so baseband noise power is «r^/4. Next we notice that white Gaussian 
noise, having same total power as quantizing noise over the band 
(— 1/2 7\ , 1/2 7\) Hz, has two sided spectral density 

N(u) = N /2 = T2'Ti watts per cycle. (35) 



* Bi is the bandwidth of gi(t),i = x, y. 
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Thus the Gaussian noise power out of the receiving filter G would be 

~- - ^ t r & i n( \ i-' 

a 2 r M 

= 12 Ti .L * (0 ~ rf/ 

+ J ^gMuMsm^Jdtj' (36) 

Since neither g x (t) nor 0„(f) are assumed to have any frequencies as 
high as w c the last two integrals above vanish. The final remark that 
completes the proof is 



/: 



9 %t)dt = o,m. 



Thus noncompandored quantizing noise behaves, at least concerning, 
its power, as zero mean white Gaussian noise, flat over the band 
(- 1/27*1, 1/27 1 ,) Hz, and total power A 2 /12. This statement is true 
with or without demodulation. 

One would like to go further and treat the baseband noise as zero 
mean Gaussian of variance given by (32). There is a justification 
for making this additional step. Recall the result of Section IV, which 
stated that if 

z = 2A, 

(2) - 2<A«) (37) 

o- 2 = varz = tV2A 2 

is a sum of independent and uniformly distributed variates A, , then 
(provided a 2 = 1/122A; < ») for all A such that A > S(A,>, 

Prob (z > A) ^ (ir)»[(p)* + (p + 2)*]P„(A). (38) 

In (38), ( P )% = 4/[(2)*<r], and P„(A) is the probability that a 
Gaussian variate of the same mean and variance as z is greater than 
A. Since P g (A) depends exponentially on p, the coefficient structure 
in (38) is not nearly as important as P g {A). We would like to argue 
(but not prove) that ignoring the coefficient in (38), that is, simply 
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assuming Gaussian behavior, is quite accurate for the baseband noise 
(30) for the error rates of practical interest. 

Thus consider eight-level, 50 percent roll off, transmission over 
our hypothetical "noncompandored 7Y' transmission facility. From 
(30) and (31) 

n h {t) = i? E «(*T. + t)K& ~ *r» - r) 

•cos [(w c - «0(* - kT, - r) + U e (kTr + r)], (39) 

where, according to an appropriate normalization, 

4r C ° S T 



-(f 



Notice that since h x « 1/i 2 , J large, the sum in (39) is bounded. 
A computer study of (39) for various values of t and t shows this bound 
to be not too sensitive (about 5 percent variations) to choices of t and t. 
Numerically we find 

l-^lsff&W (41) 

From the sum formula (26) the variance of (39) is obtained. We 
calculate 

^ - 4 12 T (42} 

Thus a peak-to-mean square ratio of the baseband noise power may 
easily be shown to be 15 dB. To obtain some insight from this value, 
consider the question of how many (N) identical independent, zero- 
mean uniform densities one would have to convolve to get a peak to 
rms value of 15 dB; the answer is N = 10. Ten uniforms generate, 
we feel, a reasonable approximation to a Gaussian curve. As a check, 
consider that our ten uniform densities each have range [—0.5, 0.5]. 
To check (not prove) the approximation on the tails we calculate 

Prob [sum ^ 4.5] = jjrj = 2.75 X 10 -7 . 

The Gaussian assumption gives 4.46 X 10" 7 . Thus we will assume 
that for error rates > 10 -7 the Gaussian assumption will yield rea- 
sonably accurate answers, not just being a bound in the sense dis- 
cussed above. 
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The above theory showing that noncompandored quantizing noise 
may be considered to be additive white Gaussian noise with zero 
mean, variance A 2 /12, of 4 kHz bandwidth has been compared with 
the experimental results of Gustafson on the performance of the VSB 
(203) data set which operates at 5400 bits per second. Fortunately 
an experimental curve is available for error rate versus signal-to- 
noise ratio without companding and this is shown in Fig. 3 along 




O -10 -20 -30 -40 

AVERAGE INPUT POWER IN dBm 

Fig. 3 — Comparison of experimental and theoretical error rates for one link 
of noncompandored transmission. Theory neglects overload distortion. Instan- 
taneous input power of +6 dBm is the onset of overload. 

with the results of present theory* (for one link of Tl.) The rise in 
the experimental curve at high input signal power results from over- 
load distortion of the quantizer which has been neglected for the 
present analysis. Overload occurs at a peak power of +6 dBm on the 
scale used in Fig. 3, and thus the peak power to average power for 
the eight-level VSB set (including pilot tones) appears to be around 
11 dB. In general the observed error rate is higher than the theoret- 



*To model the performance of the actual 203 receiver, an additional noise 
source is included at baseband, as suggested by Saltzberg 8 and shown in Fig. 2. 
The baseband S/N for this noise is chosen to be 28.08 dB. This noise alone 
would yield an error rate of 25 X 10"°. 
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ical prediction, and for low P e is even higher than the theoretical 
upper bound. Perhaps this is caused by other distortions in the system 
not considered here. 

For multiple links of transmission one should take the quantizing 
noise to have the same properties as above, but the total noise power 
;« N • A 2 /12, where N is the number of links. 

We wish to emphasize that the curve drawn in Fig. 3 does not 
represent any theoretical ideal; we have tried to understand the 
performance of an existing data set and its imperfections. Actually, 
if the data set were functioning perfectly, there is a range of input 
power where no errors would be made. We have normalized units so 
that no errors can be made if the baseband noise is less than unity 
in magnitude. Table I shows the peak value of quantizing noise cal- 
culated from equation (41) as a function of input power measured 
in the same units as in Fig. 3. An input power of —15 dB would be 
near typical operating levels. If the data set were imperfect but the 
eight-level eye were no more than 30 percent closed (but one had 
perfect timing), then one would still not make errors down to —15 
dBm. In general we see that quantizing noise is not a basic limiting 
factor on the error rate for all input power levels. 

VI. ANALYSIS WITH COMPANDING 

Equation (9) indicates that the derivative of the compressor char- 
acteristic is an important quantity. For the logarithmic curve given 
in equation (1) , 

*"(*) = i TTx \ i i 1 u. r < 43 ) 

In (1 + p) 1 + m | ;i • | 
The average of A 2 (£) now is not A 2 /12 but is 

- 5 P" V" M) ]' (1 + a»<l»l) + »'P), (44) 

where the average power P = (re 2 ). Now (\x\) cannot be less than zero 
nor more than (P) % . Hence 

*(1 + j/P) ^ <Tav ^ *[1 + K^Y (45) 

where 



IT. 



A* 
K= 12 



ML+m)] 2 . (46) 
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The lower and upper bounds in (45) indicate that for large n(P) h 
the average noise power is not a sensitive function of the probability 
density of the input signal. The knowledge of a 2 „ cannot be used here 
to obtain a strict upper bound for the probability of error as was done 
in the uncompandored case, for the "instantaneous" noise variance is 
correlated with signal values. Thus large input signals "see" bigger 
step sizes, in effect, than smaller inputs would. One concludes from this 
that for multilevel transmission the outermost levels would have the 
greatest noise associated with them. 

To make exact calculations on this matter is a difficult task, and we 
confine ourselves to some estimates of the effects. Estimates can be 
obtained by restricting attention to special sequences. Thus for an 
eight-level PAM system let an arbitrary sequence consisting only of 
the outer levels (±7) be transmitted, and compare this with another 
sequence consisting of (±5) transmitted in place of (±7). Then the 
quantizing noise will be — considering the \i term in (44) to be of prin- 
cipal importance — in the ratio 7 2 /5 2 . Thus the outer level will have, 
in this circumstance, 3 dB more noise than the next inner level. The 
contrast between these levels will be somewhat lessened in a random 
sequence using all levels, but it is clear that the 3 dB number quoted 
here provides an upper bound to the difference. 

Worst case estimates of error rate in the compandored case may be 
made by replacing A 2 / 12 in (42) by the upper bound for <j\ v given in 
(45), and finally using peak power instead of average power in (45). 

For the eight-level VSB system considered previously, operating 
on T\ facilities, this procedure yields error rates of 10~ 6 — 10~ 8 over 
one transmission link (for the interesting ranges of input power). 
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APPENDIX A 

Derivation of Quantization Error 

Equation (9) of the text relating the output sample error A(£) to 
step size A, compressor characteristic F(x), and signal amplitude 
l(t) at time of sampling is easy to derive if the chain rule is used 
to differentiate the relation 

F- l [F(x)] = x (47) 
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to obtain 

dF; v \ 



du 



u-F(i) 



dF(v) 
dv 



= 1. 



(48) 



Now clearly, if the error made when ko mp is quantized is small, that is, 
if A is small, then 



output , (/ ) , A ^ 

sample UcomiJ ^ du 



(49) 



or 



(50) 



output __ j _i_ A 
sample F (J) 

where (48) has been applied to (49) to obtain (50) . 

APPENDIX B 

Derivation of a b 

Squaring (30) and averaging over {e} gives 

•cos 2 M2 7 , + co c r] 4- gl(t - IcT, - r) 

•sin 2 [« e M\ + a> c r] + 2g x (t - kT, - r)g v (t - kT x - r) 

•sin M7\ + w c r)] cos [u € kTi + f) e r]\ (51) 



or, 



*i(t)=i(%gzi 9 i+ 9 i 



+ | (f ^)j E (A - dl) cos [2».M\ + 2a, e r] 

+ 2 £ g x g y sin [fco^T, + 2u c r]} . (52) 

All the sums in (52) may be evaluated using the Poisson sum formula 
quoted in equation (26). The first term on the right of (52) is simplest 
to handle. Since g t has no frequencies higher than 1/2 7\ , the Fourier 
transform G 2i (u>) of g 2 (t) has support contained in [— 2ir/T 1 , 2ir/T x ] } 
and further, since it is a convolution, G a <(±2ir/Ti) = 0. Thus 

| (f Q E (£ + 9l) = l (l QtAG„(0) + <?„(0)]. (53) 
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The other sums in (52) are all zero, for a typical sum is (where 6 = 
X) /(/ -nT - t) cos [nd + <p] 



= cos 



i l ¥ x + *• " yA ^ /(/ " nT ~ t) cos (/ ~ nTi ~ T) k 



(54) 

The sum formula is now directly applicable to the functions 
/(<) cos 2<a e t and j(t) sin 2w e i. The functions have Fourier transforms 
which, according to the discussion following equation (29), vanish at 
a = zk2ir/T x -k, where k is any integer, including zero. The results 
claimed follow. 
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